Temperature effects in the nuclear isoscaling 
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The properties of the nuclear isoscaling at finite temperature are investigated and the extent 
to which its parameter a holds information on the symmetry energy is examined. We show that, 
although finite temperature eff'ects invalidate the analytical formulae that relate the isoscaling pa- 
rameter a to those of the mass formula, the symmetry energy remains the main ingredient that 
dictates the behavior of a at finite temperatures, even for very different sources. This conclusion 
is not obvious as it is not true in the vanishing temperature limit, where analytical formulae are 
available. Our results also reveal that diff'erent statistical ensembles lead to essentially the same 
conclusions based on the isoscaling analysis, for the temperatures usually assumed in theoretical 
calculations in the nuclear multifragmentation process. 
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PACS numbers: 25.70.Pq, 24.60.-k 
I. INTRODUCTION 

The scaling law obeyed by the ratio of the yields 
YA.z{i) of fragments of mass and atomic nmnbers A and 
Z observed in different reactions, labeled T' and '2', i.e. 
the nuclear isoscaling P, Q 



i?2 



yA.z(2) 

Ya,z{1) 



Cexp(a7V + /3Z) , 



(1) 



opens the possibility of accessing the properties of nu- 
clear matter far from equilibrium. The parameter C in 
the above Eq. is just a normalization factor, but a and 
(3 hold valuable information on the nuclear interaction 
S 13 ■ The relevance of this property is due to the fact 
that, if the sources formed in the intermediate stages of 
the two reactions have the same temperature T , effects 
associated with the decay of the primary hot fragments 
should not distort this scaling law. Since it is difficult 
to tune the temperature of the sources experimentally, 
similar reactions are usually used so that one expects T 
to be very close in both sources. Therefore, the param- 
eters a and /3 preserve information related to the stages 
at which the fragments are created. This assumption is 
supported by theoretical calculations [1, H, H] . 

In the vanishing temperature limit and for similar 
sources, it has been shown that the symmetry energy 
coefficient of the nuclear mass formula_j_^C^yni, gives the 
main contribution to the parameter a 



T 



Ai 



Z2 
A2 



(2) 



where Ai and Zi stand for the mass and atomic numbers 
of the sources formed in the two reactions. This close 
connection between a and Csym bas extensively been ex- 
ploited in manystudies 0, 0,11, i, M, M, [13, H Q M, 

The density dependence of Cgym is of particular inter- 
est |10,i,[l3,[ii,[li,[23l- Since a is expected to retain 
the memory from the breakup configuration, it could pro- 
vide relevant information on the freeze-out stage. How- 
ever, the formal derivation of the relationship between a 
and Csym , Eq. ^ , has been carried out in the framework 
of statistical models, in which the created fragments are 
at the saturation density [1, 3- As a matter of fact, it 
has been only very recently that the thermal expansion 
of the fragments at the freeze-out stage has been consis- 
tently incorporated in the statistical models [2l|. Fur- 
thermore, although the derivation presented in Ref. [J 
clearly shows that the leading term relating a to Csym is 
indeed given by Eq. it assumed that the temperature 
of the decaying source tends to zero. Since in this tem- 
perature range the system's density must be close to the 
saturation value, it is not clear whether the conclusions 
implied by this Eq. should still hold at higher temper- 
atures, where effects associated with the multifragment 
emission are non-negligible. Nevertheless, a strict inter- 
pretation of this expression has been adopted in many 
works 0, H, [3 to explain the behavior of a at rela- 
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tively high temperatures, which led to controversial in- 
terpretations based on different scenarios [H, [3, [3, [11] . 

In this work we study the consistency of Eq. ([2]) us- 
ing the traditional Statistical Multifragmentation Model 
(SMM), in which the fragments' densities at the freeze- 
out configuration correspond to their saturation value 
[13, m, ■ In the context of this model, the expanded 
breakup volume of the system is associated with the sepa- 
ration between the fragments, rather than to the thermal 
dilatation of their volumes. Therefore, in this framework, 
the symmetry energy is strictly insensitive to the freeze- 
out volume of the system. 

We investigate whether the conclusions based on Eq. 

t, formally derived in the vanishing temperature limit 
, remain valid in the temperature range usually as- 
sumed in the nuclear multifragmentation process. Al- 
though there are no analytic formulae relating a to Cgym 
at non- vanishing temperatures, we show numerically that 
the symmetry energy rules the behavior of a at these tem- 
peratures. This conclusion remain valid also in the case 
where the other terms of the mass formula, disregarded 
in the derivation of Eq. dominate its behavior at 
low temperatures. We demonstrate that the main con- 
clusions based on this expression still hold although Eq. 
([2|) breaks down in this temperature range. More specif- 
ically, we show that, contrary to what is expected from 
Eq. 0, Csym ^ ra/4[(Zi/^i)2-(Z2M2)^], except when 
T — > 0. It should be noticed that even if former studies 

I, i, S S 1 M [m, M M M m, M [13, M, M M 

clearly show that the isoscaling parameter a and Cgym 
are strongly correlated, the extent to which the functional 
dependence given by Eq. ([2]) is a good approximation to 
the actual relationship between a and Cgym is not clear. 

The sensitivity of the isoscaling analysis to the sta- 
tistical ensemble employed in the calculation is also in- 
vestigated and we demonstrate that, in the temperature 
range relevant to nuclear multifragmentation, the micro- 
canonical (M.C.), canonical and the grand-canonical 
(G.C.) ensembles lead essentially to the same conclusions. 

The remainder of the paper is organized as follows. We 
briefly recall the main features of the SMM model in sect. 
im The results are presented in sect. IIIIl We summarize 
in sect. IIVI what we have learned. 



II. THEORETICAL FRAMEWORK 

The Helmholtz free energy F is the main ingredient of 
statistical models since, besides allowing for the calcula- 
tion of the thermodynamical properties of the system, it 
holds the physical ingredients used in the model. In the 
SMM, F is written as: 



FiT) 



C, 
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ii + xyf' ^1/3 



^trans {T) 
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-BA^z + fl.ziT)- 



where 



Ft, 



-T{M - 1) log iVf/Xl) + Tlog (go4^') (4) 
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tJ^Na^z log (5^,2^'/' 



N 



A,Z 



log(iVA,z!) 



In the above expressions, V-^/Vq = 1 + x denotes the 
ratio of the freeze-out volume to the ground state 
volume Vq of the source, Aq and Zq are its mass and 
atomic numbers, respectively. In this work, we use x = 2 
in all the calculations. The fragment multiplicity in each 
fragmentation mode is represented by Na,z, Vf = Vj^ — Vq 
is the free volume. At = ^2Trh? /mT corresponds to the 
thermal wavelength, where m is the nucleon mass. For 
A < A empirical values for the spin degeneracy factor 
dA.z sie used as well as for the binding energy B^.z- 
These light fragments, except for the alpha particle, are 
assumed to behave like point particles with no internal 
degrees of freedom. In this work, we adopt the simple 
Liquid Drop Mass (LDM) formula used in Ref. [l^ to 
calculate the binding energies of heavier nuclei: 



Ba.. 



c^,a-c.av^-c.^^+c,^ 



(5) 



where 



A-2Z 
A 



(6) 



and i = V, s corresponds to the volume and surface terms, 
respectively. For the sake of clarity in the subsequent 
calculations discussed in sect. IIIIl we have suppressed 
the pairing term. The coefficient Cc corresponds to the 
usual Coulomb term, whereas corrections associated with 
the surface diffuseness are taken into account by the 
factor proportional to Cd- Two simple versions of this 
LDM formula are used below, labeled LDMl and LDM2 
as in Ref. [l^. In the first one, fc^ = so that sur- 
face corrections to the symmetry energy are neglected 



and one has E. 



(1) 



Csym{A - 2ZflA, a 



sym 



The Coulomb correction proportional to Cd is also sup- 
pressed in the LDMl formula, i.e. Cd = 0. In the case 
of the LDM2 these terms are preserved and one would 
have siyL = [Csym - asks/A^/^] {A - 2Zf/A . For the 
sake of the comparisons below, we keep the definition 
of Csym used in the LDMl and do not include the term 
-asks/A^/\ 

The internal free energy of the fragment f\ ^ (T) is 
given by the standard SMM expression [2^ : 



fl.ziT) 



Co 



A + /3oA^/^ 



rp2 _ rp2 

C 

T2 + T2 
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, (7) 
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where /3o = 18.0 McV, = 18.0 MeV, and eq = 
16.0 MeV. This formula is used for aU nuclei with A > A. 
In the case of the alpha particles, we set /3o = in or- 
der to take into account, to some extent, the large gap 
between its ground state and the first excited state. The 
spin degeneracy factor for these nuclei is set to unity, due 
to the schematic treatment of their excited states. 

We stress that, since the parameters of the mass for- 
mula are temperature and density independent, the sym- 
metry energy of the fragments are not sensitive to the 
freeze-out density of the source. Furthermore, since in 
the version of SMM employed in this work the internal 
free energy of the fragments, Eq. ([7]), is isospin indepen- 
dent, the entropic terms of the Helmholtz free energy do 
not lead to any contribution to the symmetry energy. 

In the translational contribution to the free energy 
-Ftrans, thc factor AI — 1 (rather than AI = J2a z ^a,z) 

and Tlog^goAg^^^ arise due to the subtraction of the 

center of mass motion. 

These ingredients are used in all the different ensembles 
employed in this work, i.e., micro-canonical, canonical, 
and grand-canonical. For practical reasons, as it allows 
for extremely fast calculations, we used the McGill ver- 
sion of the canonical ensemble of SMM [l^. Since we 
adopt the same ingredients, it is equivalent to the tradi- 
tional canonical Monte Carlo version of SMM [2^. We 
refer the reader to Ref. [l3| for details on the grand- 
canonical ensemble calculations. 



III. RESULTS AND DISCUSSION 

As has been demonstrated several times [1, 0; @1 ; the 
isoscaling parameter a is associated with the baryon 
chemical potentials /^^\ i = 1,2, of the sources Si and 
S2 through 



(2) 

Ms 



T 



(8) 



This relationship is obtained in a scenario in which the 
fragments' densities correspond to the saturation value 
and, therefore, simple analytical expressions for the av- 
erage yields can be obtained [lB|. Otherwise, the free vol- 
ume dependence on the fragments' multiplicities [2l| lead 
to highly non-linear terms in the Helmholtz free-energy. 
As a consequence, the average yields are no longer given 
by the simple traditional grand-canonical formulae. Since 
in this work the fragments' densities are fixed and corre- 
spond to their ground state values, these difficulties can 
be disregarded. One of the underlying assumptions made 
in the derivation of Eq. ([5]) is that the two sources have 
the same temperature. 

The chemical potentials are strongly influenced by the 
LDM employed in the calculation as well as by T and the 
composition of the source. In the vanishing temperature 
limit, it has been demonstrated [J for the LDMl that 



0) = -a„ + - 



.1/3 



-a.- 



' A 



4/3 



-a 



sym 



1 - 
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where — Cc[l — 1/(1 + H is easy to extend this 

formula to the LDM2, also in the T ^ limit, and from 
Eq. IS]) one obtains: 
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In the above expression, we have introduced the short- 
hand notation A(Z/^)2 ee \{ZllA\ - ZljAl). 

It should be noted that in the statistical models the 
properties of the sources are taken for granted, as the 
dynamics which leads to the configuration at the freeze- 
out stage is beyond the scope of the approach. Therefore, 
we explicitly assume that the sources S\ and 6*2 have the 
same temperature T . The validity of this assumption 
must be investigated on the basis of dynamical models. 
For the present purpose, which consists in verifying the 
consistency of Eq. ^ with the statistical treatment upon 
which this expression is based, the pair of sources (Sx^S'i) 
are conveniently selected according to their Zj A ratio, as 
well as their sizes. 

If the mass and atomic numbers of the sources are 
similar, the symmetry energy dominates the sum and, 
by neglecting factors involving a^, ajts-, ^ well as thc 
Coulomb terms, it leads to Eq. (I2|), usually employed in 
the literature. From this, one may obtain the apparent 
symmetry energy coefficient: 



Unc 
sym 



Ta/A(Z/A)2 



(11) 



which we label as "Unc" to emphasize that the correc- 
tions associated with the remaining terms of the mass 
formula have been discarded. 

If the sources S\ and S2 are not similar enough, i.e. 
the differences of their mass and atomic numbers are not 
small, the corrections neglected in the above expression 
introduce a fixed shift to Cgym for all temperatures, whose 
value depends on the sources considered: 



C£" = {Ta^ C[. 
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FIG. 1: (Color online) Apparent symmetry energy coefficient 
C^yl^'^°' for different systems (with A(Z/yl)^ = 0.147) cal- 
culated with the G.C. ensemble for the LDMl formula. The 
dashed lines correspond to a^fc„ = 22.39 MeV. For details, 
see the text. 



and the label "Cor" denotes this fact. 

Since the above expressions, Eqs. PT|) - P^ . rely on the 
baryon chemical potential they must be accurate only in 
the T — > limit. Therefore, it is not clear whether the 
conclusions implied by any of them remain valid at finite 
temperatures. The extent to which a remains related to 
Cgym at these temperatures is not a priori obvious. It 
should be observed that studies based on the statistical 
|Jii,Bj9|, [13, [m, E, [H, Q m, H [3 and dynamical 
la,[l3j[20| models have only shown that a is correlated to 



FIG. 2: (Color online) Same as Fig.[T]for the LDM2 formula. 
The dashed lines correspond to a«A;„ = 27.80 MeV. For de- 
tails, see the text. 



Cgym but the accuracy of Eq. ^ remains unknown. Since 
there are no analytical expressions for at finite tem- 
peratures, for configurations in which many fragments 
are present, this issue is numerically investigated below 
in the framework of the SMM. 

In this way, we reconstruct Ci^m^*^°' through Eqs. 
(fTT|) - (|12[) and check whether it is consistent with the 
fixed parameters used in the LDM adopted in the model. 
The parameter a is calculated from fits based on Eq. 
(P) by the yields provided by the G.C. ensemble (except 
where stated otherwise) in all the calculations presented 
below. If either of the above Eqs. is valid at finite tem- 
peratures, by construction, one should obtain a constant 
value Csyni = ciiifc-u for all T . Indeed, since the low density 
values used in the present SMM models is associated with 
the spacing among the fragments, rather than with their 
thermal volume expansion, the reconstructed C[l^m^^°' 
should be constant and it should be given by a^ky for 
any temperature and selected sources, as these parame- 
ters are fixed at the saturation density value for all nuclei 
[isl. Therefore, deviations from the constant value 
(Xd kfj SI gnal the breakdown of the exact relationship be- 
tween a and Cgym implied by Eqs. (jlip - (|12p and usually 



Unc/Cor 



as a func- 



assumed in isoscaling analyses. 

Figure [T] displays the apparent C. 
tion of the temperature, obtained with the expressions 
above, for the LDMl {k, = 0) and the ("2Sn,i24Sn), 
(i^8Re,i24sjj)^ (56Mn,62i^n), and (^^^Rc^^Mn) pairs of 
sources {Si, 82)- One should notice that A{Z/A)^ is fixed 
for these sources, so that the main contribution associ- 
ated with the symmetry energy is the same in all the 
cases. The results corresponding to Eq. (|lip are dis- 
played in the upper panel of this picture and are very 
similar for the ("2Sn,i24Sn) and (5<'Mn,62Mn) sources. 
As expected, approaches a^ky — 22.39 MeV at low 

temperatures, but important deviations are observed in 
this limit for the (i6*Re,i24Sn) and (i68Re,S2]vin) sources. 
This is because the contributions associated with as and 
the Coulomb terms in Eq. (|10p cannot be disregarded in 
these cases. The results displayed at the bottom panel 
of Fig. [H C^°^, show that the expected value of Cgym 
is obtained in the low temperature limit when the corre- 
sponding shifts, Eq. (|12p . are taken into account. 

One striking feature of the upper panel of Fig.[T]is that, 
in spite of the large differences observed at low tempera- 
tures between the (i^^Re,i2''Sn) or {^^^Re,^^Mn) sources 
and the others, the four curves merge for T > 5.0 MeV. 
This means that Eqs. (fTTI) and (fT^ break down at large 
temperatures and a can no longer be approximated by 
Eq. ([2]) or (fTO]) . However, the fact that the curves given 
by Eq. ([TT]) agree for T > 5.0 MeV implies that, although 
this Eq. is not valid in this regime, a is indeed strongly 
correlated with Csym, despite the fact that the analyti- 
cal temperature dependence of a is unknown. Further- 
more, it also shows that the Coulomb and surface terms 
play a negligible role for T > 5.0 MeV, although they 
cannot be neglected at smaller temperatures. Since the 
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FIG. 3: (Color online) Isoscaling parameter a calculated with 
the LDMl formula. For details, see the text. 
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nuclear multifragmentation process is expected to take 
place from this temperature value on, these results sug- 
gest that a remains a valuable observable in these studies. 

In order to check whether these conclusions are not bi- 
ased by the simple form of the LDMl formula, we also 
performed calculations using the LDM2 model. The re- 
sults are shown in Fig. [21 One sees that the differences 
are much more pronounced in this case at low tempera- 
tures but the qualitative conclusions above remain true. 
For T > 5.0 MeV, is the same for all the systems, 

in spite of the very large differences in the vanishing tem- 
perature limit. This suggests that our findings will not 
be affected by the use of more involved mass formulae. 

The above results show that the relationship between 
a and Csym is much more complex than is suggested by 
Eqs. ([2]) or pO)) . The constraints to which the chemical 
potentials are submitted at large temperatures, where 
many fragments are present, lead to important changes 
to Eq. ([5]), on which these Eqs. rely. In the particular 
case of the SMM, these results suggest that 



X ao (1 + aiT -f- aaT^ -f agT^) 



(13) 
(T > 5.0 MeV) 



For the LDMl, one has oq = 1.3338, ai = -0.131842 
MeV-\ as = 0.0156409 MeV^^, and 03 = -0.000531862 
MeV-3, whereas oq = 2.6431, ai -0.292421 MeV-\ 
02 = 0.0338321 MeV-2, and 03 = -0.00126698 MeV"^ 
in the case of the LDM2. The values of Cgym recon- 
structed from this expression are depicted by the squares 
in the upper panel of Figs. [1] and [2] The results show 
that Csym is fairly constant and that its value is con- 
sistent with the parameters k^a^ employed in the LDM 
models. However, the form of the above expression has 
been chosen for illustration purposes and the actual re- 
lationship between a and Cgy-m remains unknown. To 
obtain a closed expression, one must be able to derive 



FIG. 4: (Color online) Uncorrected apparent symmetry en- 
ergy coefficient calculated from the LDM2 formula for 
£0=4 MeV (upper panel). The bottom panel shows a com- 
parison between the results obtained with eo = 4 MeV and 
eo = 16 MeV. For details, see the text. 



analytical expressions for the chemical potentials associ- 
ated with a many fragment system at finite temperature, 
which has not been achieved yet. 

To sec that the symmetry energy is, in fact, the essen- 
tial ingredient that dictates the behavior of the isoscal- 
ing parameter a at high temperatures, we consider the 
(i68Re,ii2Sn), (ii2Sn,56Mn), and (i68Re,56Mn) sources, 
for which A{Z/A)'^ = 0. Equation (dU]) predicts that 
a due to contributions from the Coulomb and surface 
terms, whose relevance increases as the differences be- 
tween the mass and/or the atomic numbers of the sources 
become larger. The results displayed in Fig. [3] show that 
this is true only for T — > 0. The effect of these terms van- 
ishes at high temperatures and a — > for T > 5.0 MeV, 
i.e. one sees once more that the Coulomb and surface 
terms have very little influence on the behavior of a at 
high temperatures. The symmetry energy is the key in- 
gredient in this temperature domain, despite the formu- 
lae derived at T ^ 0. 

Since there are no analytical expressions for a at fi- 
nite temperatures, we investigate the influence of the in- 
ternal excitation terms of the free energy by increasing 
the level density of the bulk term of Eq. ([7]) and we set 
£0=4 MeV. One should notice that this leads to an 
important increase in the magnitude of f*{T). For in- 
stance, for A = 100 and T = 6.0 MeV, the standard 
value is f*{T) = -319.4 MeV whereas for eo = 4 MeV 
one has /*(T) = -994.4 MeV. The apparent C^^^^ is 
shown in the upper panel in Fig. [5] for eo = 4 MeV and 
the LDM2 formula. One sees that the qualitative fea- 
tures observed previously remain the same and that all 
the curves merge for T > 5.0 MeV. Furthermore, the re- 
sults displayed in the bottom panel of this picture show 
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FIG. 5: (Color online) Isoscaling parameter a calculated with 
the M.C, G.C and canonical ensembles for the LDMl and 
LDM2 formulae. For details, see the text. 



that, in spite of the large differences in /*(T), the effect 
on the apparent Cj^"^ is very small. This reveals a fair 
insensitivity of a to the entropic terms of the internal 
free energy of the nuclei. We have checked that similar 
results are obtained if one reduces Po in Eq. which 
gives the surface contribution to f*{T). 

Finally, we check whether our conclusions depend on 
the statistical ensemble employed in the calculations. 
Figure [5] shows a comparison between the values of a 
obtained with the M.C, canonical, and G.C. ensembles, 
using the two LDM formulae discussed in this work. One 
observes very small differences at large temperatures, i.e. 
for T > 4 — 5 MeV, which is the temperature domain in 
which these models are usually employed in the study 
of the nuclear multifragmentation. The discrepancies at 
low temperatures between the canonical and the G.C. en- 
sembles are due to finite size effects, since these systems 
are far from the thermodynamical limit psj . The devia- 
tions are more important in the case of the M.C. calcu- 
lations at low temperatures where the excitation energy 
of the system is very small and the strict energy conser- 
vation constraint plays an important role in determin- 



ing the most important partitions. These results suggest 
that our conclusions should not be affected by the use of 
different ensembles in the temperature domain which is 
relevant to the multifragment emission. 



IV. CONCLUDING REMARKS 

We have demonstrated that the isoscaling parameter 
a is not sensitive to the Coulomb and the surface terms 
of the nuclear binding energy in the temperature domain 
in which the multifragment emission is expected to take 
place, although they are very important in the vanishing 
temperature limit. The symmetry energy dictates the 
behavior of a for T > 5.0 MeV, despite the fact that the 
simple analytical formulae that relate a to Cgym, Eqs- © 
or pO|) . are not valid in this temperature domain. The 
weak dependence of a on the entropic terms of the in- 
ternal excitation energy of the fragments suggests that, 
at high temperatures, a is essentially governed by the 
symmetry energy and T, although its actual functional 
dependence is rather complex in this regime. More pre- 
cisely, we found that a is given by an unknown func- 
tion / which, for T > 5.0 MeV, is not sensitive to the 
Coulomb and surface terms, i.e. a = /(Csym, (-^iMi)^ — 
{Z2 / A2)'^ , T) . Furthermore, we also found that a 
4CsynJ(^iMi)2 - (Z2M2)2]/T, except for T ^ 0. De- 
spite this fact, since a = /(Csym, (Zi/Ai)'^ -{Z2/A2)'^ ,T) 
in the temperature domain at which the multifragment 
emission is expected to occur, it suggests that a is indeed 
a good probe for the symmetry energy coefficient but in- 
terpretations based on the isoscaling analysis should be 
taken with care. We also found that, for T > 4 or 5 MeV, 
the isoscaling analysis is not sensitive to the statistical 
ensemble employed. 
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